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Abstract 

^ . We construct the equations of motion for the coupled dynamics of order parameter and concen- 

, tration for the nematic phase of driven particles on a solid surface, and show that they imply (i) 

giant number fluctuations, with a standard deviation proportional to the mean and (ii) long-time 
' tails ^ t^"^/^ in the autocorrelation of the particle velocities in d dimensions despite the absence 

o , 

, of a hydrodynamic velocity field. Our predictions can be tested in experiments on aggregates of 

amoeboid cells as well as on layers of agitated granular matter. 
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Recent experiments [1, 2] offer evidence for long-ranged orientational order in collections 
of driven particles. The work of [1] is concerned with aggregates of hving cells confined to 
a solid surface, whereas that of [2] studies vertically vibrated layers of granular rods. Two 
kinds of uniaxial orientational order are discussed in [1], and appear to arise in [2] as well: 
(i) nematic or apolar, where the axes of the elongated cells are oriented, say along a mean 
direction n, but n and - n are equivalent; and (ii) vectorial or polar order, where n and - n 
are not equivalent. On symmetry grounds [3], a vector-ordered nonequilibrium steady state 
should have a nonzero macroscopic drift velocity Vq, while for true nematic order Vq = 0. 

We shall be concerned here with strictly nematic order. In [2] a nematic phase is men- 
tioned but not discussed much. In [1], there is particularly convincing evidence for a nematic, 
including images of classic strength -1/2 disclinations, from studies on melanocytes (cells 
which distribute pigment in the skin). In fact in this class of cells the individual particles 
themselves are head-tail symmetric, or apolar, so that the drift velocity is zero for each par- 
ticle. The cells involved are tens to hundreds of microns in length, and the rods in [2] are 6 
mm long, so that thermal Brownian motion plays no role in the dynamics of these particles. 
What distinguishes the particles in these experiments from conventional nematics is their 
active or driven character: they absorb energy from their surroundings (from nutrient in 
one case, shaking in the other) and transform it into motion in the plane. The dynamical 
behaviour of melanocytes, for example, involves rhythmic movements of the cell body as 
well as two long projections called dendrites which appear to be responsible for the inter-cell 
interaction and the level of fluctuations [1] . In granular rods the interaction and noise come 
from coUisions. Ordered phases of particles on an air-table [4] would also obey our model, 
except for possible complications arising from the periodic array of air-holes. 

As has been noted earlier [5, 6] the organization, dynamics, and fluctuation properties of 
such active, energy-dissipating objects cannot be described by equilibrium statistical physics. 
Continuum models appropriate for the dynamics of polar order wherein the particles are on 
average in motion were presented in [6]; the extension of those models to include momentum- 
conserving hydrodynamic flow and both polar and apolar order was discussed in [7]. The 
cells in [1] and the rods in [2], however, move on a solid surface which acts as a momentum 
sink for the velocity even at zero wavenumber. 

In this paper we derive the most general universal equations of motion, valid for all 
systems of active nematogenic particles without total momentum conservation [8], for the 
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apolar order parameter and concentration fields in such systems, retaining all relevant terms 
allowed by symmetry. Our aim is to look for universal properties, independent of micro- 
scopic details, that distinguish these active nematics from their counterparts at thermal 
equilibrium. 

Two predictions stand out as most striking: (i) The standard deviation in the number 
N of particles is enormous, scaling as in the entire nematic phase for two dimensional 
systems. By contrast, in any thermal equilibrium system with a finite compressibility, which 
in practice means any equilibrium system not at a continuous phase transition, number 
fluctuations are far smaller, scaling as \/N as ^ oo. (ii) The autocorrelation of the 
velocity of a tagged particle decays with time t as t^''/^ for space dimension d > 2. It 
should be recalled that precisely this form of long-time tail [9] occurs in bulk fluids at 
thermal equilibrium, but as a result of advection of particles by thermal fluctuations in the 
hydrodynamic velocity field, which decays slowly because of momentum conservation. Here 
the total momentum of the particles is not conserved, since it is damped by friction with the 
substrate [10]. A collection of particles at thermal equilibrium on a substrate would show 
velocity autocorrelations decaying much more rapidly [11]; specifically as It should 

be possible to test these predictions by studying large-scale real-space images and videos of 
the behavior of systems such as those in [1] as well as two-dimensional nematic phases in 
agitated layers of granular particles [2] . 

Our predictions for nematically ordered phases of powders stand in dramatic contrast to 
experimental results obtained on disordered, monodisperse powders (ping-pong balls) [4], 
which showed that the latter systems could be described surprisingly completely by equi- 
librium statistical mechanics, albeit with a very high Active temperature [12]. The findings 
we've summarized above clearly show that nematic phases of rod-like powders cannot be 
described by equilibrium statistical mechanics at all, even at the crudest level of scaling [13] . 

We turn next to the construction of our equations of motion. 

Since our nematogens are spread on a solid substrate, there is only one conserved quantity 
of relevance, namely, the number of particles. A description valid at sufficiently long length- 
and time-scales [14] needs then to include only the concentration field c(r, t) of particles and 
the orientational order parameter, described by the traceless symmetric tensor field Q with 
components Qij{r,t) at point r and time t. For uniaxial nematics, to which we restrict our 
attention here, we can write Q = [nn — {l/d)\]S, where the unit vector n is the director 



3 



field, I is the unit tensor, and the conventional scalar order parameter S [15] measures the 
degree of order. It is however convenient to introduce and later eliminate a fast variable, 
the velocity field v(r, t) of the particles. Number conservation tells us that 

where the number current j = cv. Newton's 2nd law for the local momentum density mj, 
where m is the mass of a particle, gives 

= -rv - V • (7 + fj? = -r V - WoVc - wiV ■ (Qc) + fn (2) 

up to bilinear order in the fields. In (2), the term in F represents friction with the substrate, 
the stress tensor a contains the effect of interparticle interactions, the Wq term can be thought 
of as osmotic pressure, is a random, nonconserving, spatiotemporally white Gaussian 
noise, not in general of thermal origin, and the wi term, unique to this driven system, says 
that inhomogeneities in the nematic order drive mass motion. 

To see where this comes from, note first that, by definition, inhomogeneities in a produce 
local acceleration. What is novel here is that there is a contribution to the stress tensor 
proportional to the nematic order parameter. This is allowed by symmetry here, and can be 
obtained by noting (see [7, 16]) that each active particle is a force dipole. Such a term is not 
allowed in an equilibrium nematic since the stress tensor in equilibrium must be derived from 
functional derivatives of the free energy. Since that free energy must be rotation invariant, 
it can only depend on derivatives of the nematic director n (in fact, this argument leads to 
the familiar Frank free energy [15] for equilibrium nematics). Functional derivatives of this 
then lead to a stress a which depends on derivatives of n. Since wc arc considering non- 
equilibrium nematics, however, this constraint that a be derived from a functional derivative 
of a rotation-invariant free energy does not apply. Hence, the only constraint on allowed 
terms is rotation invariance of the equations of motion, which is respected by a term in a 
proportional to Q [17], since both objects are symmetric second rank tensors. 

We now restrict our attention only to small fluctuations about the uniaxial nematic phase, 
aligned along the z axis, and denote vector components in directions normal to z by the 
subscript ±. We assume the nematic order is well developed, so that S can be treated 
as constant, and consider spatial variations c{r,t) — cq + 5c{r,t), where cq is the mean 
concentration, and n(r, t) = z + 5n(r, t). Since n is a unit vector, we need consider only the 
transverse components ^n^. 



To close our equations, we need the equation of motion for n [7], which is identical in 
form to that for an equilibrium nematic and reads, to linear order: 

dtSn^ = X+d,vi_ + X-\/±v, + KiV±{V± ■ 5n^) + XsVi^n^ + KsS^Sn^ + (3) 

Here A± = (A ± l)/2 ; A is the "flow alignment parameter", familiar from equilibrium 
nematics, which affects the response of the nematic director to shear, and the noise is 
delta-correlated in space and time: 

{f±iir, t)hjiv', t')) = 5ij Aj'{r - r')5{t - t') (4) 

We now use the equation of motion (2) to eliminate v. For hydrodynamically slow processes, 
dtV <^ Fv, and so we can neglect the dtV term in (2) and immediately solve for v, finding, 
to linear order, 

V = -a ( d^5n± + (V_l • Sn±) z) - 7iV_lc - 72 9^ cz + ^ (5) 

where a and 7^ are related to Co, S and the coefficients in (2). Inserting (5) into (1) gives, 
to linear order in fluctuations, the concentration equation of motion 

= {D,dl + D^Vl) 8c + 2co a d,{V 1. ■ 5n^) + V • (6) 

where the various coefficients in (6) and the noise fc are derived from those in (2) and (5). 
The leading nonlinearities that arise in (6) are of the form V_l • {dzSn^Sc) and similar terms, 
as well as (V_L • (5nj_)^ and similar terms. These nonhnearities as well as the couphngs to 
the director in (6) are forbidden at thermal equilibrium. 

The director dynamics obtained by inserting (5) into (3) is given to linear order by 

^ = {K,dl + K^Vl + Vx- ) + D,^d,V^5c + (7) 

where, again, the phenomenological parameters and noise source fj_ are derived from corre- 
sponding terms in (2) to (5). The leading nonlinearities which will appear in (7) are of the 
form 925nx(V_L • 5nj^), 5Ta.x^dz{^ _l ■ 5nj^) and similar terms. These nonhnearities too are 
completely nonequilibrium in origin. 

The concentration field as well as the nonlinearities mentioned above are absent in the 
treatment of Gruler et. al. [1]. The effect of the nonlinearities will be dealt with elsewhere 
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[18]. The inclusion of the concentration field leads to our most striking nonequilibrium effect, 
namely, the giant number fluctuations we predict below. 

The linearized equations of motion are easily analyzed for their mode structure, pos- 
sible instabilities, and fluctuation statistics. This is most conveniently done by Fourier- 
transforming in space and time, i.e., considering modes that go as exp(iq.r — icut). Doing 
this enables us to rewrite (6) and (7) as 



iuj + Dc (q) 6c (q, a;) + 2coaqz q± Stil (q, uj) =iq- fc (q, uj) 



(8) 



Dcnqzq±Sc (q, uj) + -iu Kl (q) q^ Sul (q, uj) = /l (q, uj) 



(9) 



where 



—iuj + Kt (q) q"^ Snx = fr 



(10) 



/l (q, uj) = q±-f±,fT^f±- q± (q± • f±) 



(11) 



5ni(q,a;) = q_L (q_L • Sn±{q,uj)) 



(12) 



and 



Sut = Sn± — Sul 



(13) 



are, respectively, the components of fj, and Sn± along and transverse to q_L [19] and we've 
deflned the direction dependent diffusion constants 



D, {q) = D,ql + D^q[ 



(14) 



KL(q) = K,qi + (K^ + K'L)q' 



(15) 



KTicD=K,ql + K^q[ 



(16) 



The eigenfrequencies obtained from these equations are all diffusive. There are d — 2 trans- 
verse diffusive modes coming from the equation, all with identical imaginary frequencies 



ujt — —iKt (q) q^ 



(17) 



and two coupled concentration - longitudinal modes with 



a;± = -ir±(q)g2 (18) 

with 

' ^ ' (19) 



r.(q)-^ 



D, (q) + Kl (q) ± ^J{D, (q) - Kl (q))' + 8coaD,nqlq 



2 ;52 
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We shall assume linear stability, which is assured if the Ki% and Di% are positive and the 
combination c^aDcn is not too large. 

With these eigenfrequencies in hand, it is straightforward to solve (8) -(10) for 5c, ^n^, 
and ^n^^ obtaining 

5c = (2co«g±g.) /l . . 

(u; + ir+ (q) + (q) g2) ^ ^ 

(a; + ir+ (q) + (q) ^ ) 



Snr = . ^ ^ 2 (22) 

-zo; + i^T (q) q 

Autocorrelating these, using the known autocorrelations Eqn. (4) for the random forces gives 

A. (o;^ + /^g(q)g^)An .... 

/ = (c.^ + n(q),^)(c.^ + r^_(q),^) 

/ / u 2«coag^g± (cj + iL'c (q) g^) A„ 

(q, c.) <5c (-q, -.)) = ^ ^4) + ^4) (25) 

Snr^ (q, c.) (-q, -.)) ^ ^^f^f^i^^^^,^ (26) 

with {SnTiSriL) — {SriTiSc) — for all i, q and cu. 

In writing eqns. (23 - 26), we have neglected all of the contributions coming from the 
noise fc in the concentration equation of motion (5), since these are all negligible, as q — > 0, 
relative to those coming from the noise in the director equation of motion (3). 
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Equations (23 - 26) can now be Fourier-transformed back to get real space-real time 
correlations functions. One quantity of particular interest are the equal-time, spatially 
fourier-transformed correlations of the concentration: 



■oo 



where 



F (q) = 2 (De (q) + Kl (q)) [d, (q) Kl (q) - 2coaD,nqlql] (28) 



The important thing about this result is that the {Sc)'^ fluctuations are enormous as 
g — > 0, diverging as ^. This would be equivalent, in an equilibrium system, to having a 
compressibility which diverged 0. To restate this in real space, it is as if the 

compressibility of a system of linear extent L diverged according to x oc as L — > oo. 
Independent of any connection to a response quantity like the compressibility, the rms 
number fluctuations \/SN^ in a volume V scale as ^S{q —>■ 0)V. For active nematics this 
implies that 

VSN^ (X Vm^ (X L^^^ (X N^+^ (29) 

where in the last step we've used the fact that N ocV (x L"^ and solved for L as a function 
of A^. For d = 2 this recovers our earlier statement that ^JJdPP) oc N. 

Lastly, Eqn. (5) for the velocity implies that the autocorrelation of a tagged particle will 
be controlled mainly by the effect of soft director fluctuations: a particle at a point in the 
nematic moves with a speed proportional to the local curvature of n^. This means the 
autocorrelation of the velocity v{t) of a particle with trajectory R(i) will be, schematically, 

(v(0) • v(i)) ~ {V5n^{R{0),0)VSn^{-R{t),t)), (30) 

Fourier transforming the n correlation functions (25) and (26) back to real space and time, 
and using the fact that, for a diffusing particle, ^|R.(^) — R(0)|^^ oc t, implies that the 
correlation function (30) decays as i"''/^, as claimed earlier [20]. 

We close with some comments on the experiments in [2] . Although the nematic "phase" 
in their images appears rather polydomain, it is likely that their system does display true 
nematic order in some parameter range. Perhaps the polydomain structures coarsen at long 
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times to a true nematic. Such a driven nematic would be the ideal place to test the predic- 
tions of this paper. If, for some range of parameters, they are able to find macroscopically 
aligned regions of tilted rods — which would constitute vectorial order and therefore, as 
explained in [2], move coherently in one direction — these would be examples of the moving 
XY-model of [5, 6]. Indeed, the coupled equations for tilt and concentration in [21], which 
analyzes the experiments of [2], are closely related to those of [6]. Lastly, the "vertical 
rods" phase, from the images, displays a good deal of order; one wonders whether it is a 
nonequilibrium realization of a crystal or hexatic. 
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